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Laser  cooling  of  bulk  matter  uses  thermally  assisted  fluorescence  to  convert  heat  into  light  and  can 
be  interpreted  as  an  optically  pumped  laser  running  in  reverse.  Optical  pumping  in  such  devices 
drives  the  level  populations  out  of  equilibrium.  Nonthermal  radiative  energy  transfers  are  thereby 
central  to  the  operation  of  both  lasers  and  luminescent  coolers.  A  thermodynamic  treatment  of  their 
limiting  efficiencies  requires  a  careful  development  of  the  entropy  and  effective  temperatures  of 
radiation,  valid  for  the  entire  range  of  light  from  the  blackbody  to  the  ideal  laser  limiting  cases.  In 
particular,  the  distinct  meaning  and  utility  of  the  brightness  and  flux  temperatures  should  be  borne 
in  mind.  Numerical  examples  help  illustrate  these  concepts  at  a  level  suitable  for  undergraduate 
physics  majors.  ©  2005  American  Association  of  Physics  Teachers. 
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I.  INTRODUCTION 

Energy  transfers  between  thermodynamic  systems  are  of¬ 
ten  labeled  as  either  work  or  heat.1  However,  it  is  generally 
recognized  that  the  exchange  of  electromagnetic  radiation 
cannot  be  neatly  categorized  as  one  or  the  other.  Exceptions 
are  blackbody  radiation,  which  is  one  of  the  three  traditional 
forms  of  heat  transfer,  and  ideal  laser  radiation  which  can  be 
characterized  as  work  by  the  Caratheodory  definition.2  Some 
authors  explicitly  distinguish  radiation  as  a  third  category  of 
energy  transfer.3  However  there  is  nothing  special  about  ra¬ 
diation  in  this  regard;  any  irreversible  transfer  of  energy  be¬ 
tween  two  systems  generally  involves  a  mixture  of  heat  and 
work.  For  example,  when  a  block  slides  over  a  rough  table, 
there  is  both  mechanical  work  involved  in  the  bending  of 
asperities  on  the  surfaces  of  the  block  and  the  table,  and  heat 
transfer  between  portions  of  the  block  and  table  at  different 
temperatures.4 

The  primary  reason  for  delineating  heat  from  work  is  to 
introduce  the  second  law  of  thermodynamics  and  arrive  at 
the  concept  of  entropy.  The  limiting  efficiency  of  an  engine 
or  refrigerator  depends  on  the  entropies  of  the  input  and  out¬ 
put  energy  fluxes.  This  paper  explicitly  treats  photon  engines 
and  optical  coolers,  thereby  requiring  a  robust  understanding 
of  the  thermodynamic  terms  that  apply  to  radiation. 

II.  BASIC  CONCEPTS  OF  LASER  COOLING 
OF  BULK  MATTER 

Before  jumping  into  the  abstractions  of  radiation  entropy 
and  temperature,  consider  the  principles  of  operation  of  an 
optical  engine  or  refrigerator.  Figure  1  is  a  simple  block  dia¬ 
gram  of  a  laser  cooler,  presented  in  the  same  style  that  phys¬ 
ics  textbooks  typically  use  to  introduce  refrigerators.  The 
left-hand  arrow  denotes  an  external  source  of  energy  Ein.  In 
a  conventional  refrigerator  this  external  source  does  electri¬ 
cal  work,  but  in  the  present  context  it  is  a  directed,  coherent 
(and  possibly  polarized,  particularly  if  the  active  material  is 
axial)  narrowband  laser  beam. 

The  refrigerator  consists  of  an  optically  active  medium 
that  resonantly  absorbs  most,  if  not  all,  of  the  input  laser 
radiation.  (The  absorption  length  may  be  increased  either  by 
shaping  the  material  into  the  form  of  an  optical  fiber  or  by 
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placing  it  in  an  optical  cavity  tuned  to  the  pump  wavelength.) 
Subsequently  the  medium  relaxes  by  spontaneous  emission 
at  a  variety  of  different  possible  energy-level  transitions,  so 
that  the  output  fluorescence  (which  carries  away  energy  Ea ut) 
is  multidirectional  (probably  isotropic),  incoherent,  unpolar¬ 
ized  (or  at  most,  partially  polarized),  and  broadened  in  band¬ 
width.  (This  fluorescence  must  be  absorbed  by  some  external 
heat  sink  which  is  thermally  insulated  from  the  refrigerator. 
Care  must  be  taken  to  minimize  reabsorption  of  the  exiting 
fluorescence,  for  example  by  employing  a  fiber  sample  ge¬ 
ometry  that  enables  the  fluorescence  to  readily  escape  out  its 
sides.)  In  ordinary  cases  of  laser  absorption,  Eout  is  less  than 
Eia,  and  the  difference  appears  as  an  increase  in  the  internal 
energy  of  the  material,  thereby  enabling  lasers  to  cut,  weld, 
or  ablate  substances.  However,  the  crucial  and  surprising 
property  of  certain  materials  irradiated  with  lasers  of  just  the 
right  color  is  that  they  emit  light  of  a  higher  frequency  (and 
hence  photon  energy)  than  they  absorb;  this  property  is 
known  as  anti-Stokes  or  thermally  assisted  fluorescence.  If 
anti-Stokes  emission  predominates  over  resonant  and  Stokes 
processes,  then  heat  Q  is  withdrawn  on  average  from  the 
medium  in  each  excitation-relaxation  cycle. 

The  simplest  possible  level  scheme  of  such  an  absorber  is 
depicted  in  Fig.  2.  The  material  has  only  three  energy  levels, 
where  level  1  is  the  ground  state.  The  pump  laser  frequency 
vp  is  tuned  to  resonance  with  the  1-2  transition,  so  that 
Em=hvp  (times  the  number  of  photons  absorbed),  where  h 
is  Planck’s  constant.  Suppose  that  vp  is  large  enough  that 
nonradiative  relaxation  from  the  excited  states  to  the  ground 
state  is  much  less  likely  than  radiative  decay.  In  other  words, 
the  fluorescence  quantum  efficiency  is  essentially  unity.  Ac¬ 
cording  to  the  energy  gap  law,5  this  requirement  is  satisfied  if 
vp>v acceptor >  where  Acceptor  is  the  frequency  of  any  other 
modes  of  the  substance  to  which  the  excited  states  can 
couple.  For  a  solid  material,  h  vacceptor=kT  (where  k  is  Boltz¬ 
mann’s  constant)  could  represent  phonon  energies,  where  kT 
is  a  typical  thermal  energy  at  the  refrigerator’s  operating 
temperature  T.  Alternatively,  the  energy  acceptors  might  be 
localized  vibrational  or  rotational  modes.  If  we  assume  that 
the  radiative  time  constant  for  relaxation  from  the  excited  to 
the  ground  states  is  relatively  long,  as  is  typical  of  fluores¬ 
cent  materials,  and  that  the  spacing  between  levels  2  and  3  is 
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Fig.  1.  Block  diagram  of  the  input  and  output  energy  flows  for  a  laser 
cooler. 


Fig.  3.  Thermodynamic  diagram  analogous  to  Fig.  1  describing  a  laser 
pumped  by  an  external  light  source. 


small  (compared  to  kT),  so  that  rapid  nonradiative  ex¬ 
changes  occur  between  them,  then  a  thermal  (Boltzmann) 
distribution  of  population  between  these  two  excited  levels 
will  be  established  prior  to  relaxation  to  the  ground  state.  In 
that  case,  some  fraction  of  the  active  species  (viz.  ions,  at¬ 
oms,  molecules,  or  electrons)  will  be  thermally  promoted 
into  level  3.  (The  fraction  depends  on  the  relative  degenera¬ 
cies  and  energy  spacing  of  levels  2  and  3.)  Consequently,  the 
average  fluorescence  photon  energy  liVt  will  be  larger  than 
the  laser  pump  photon  energy  h  vp ,  with  the  differency  rep¬ 
resenting  thermal  energy  withdrawn  from  the  system,  Q 
=  h(Vf—  vp)  per  cycle. 

This  cooling  energy  Q  may  be  used  to  chill  an  external 
load.  Typical  proposed  applications6  of  optical  refrigerators 
are  to  cool  infrared  detectors  in  satellites  (because  competing 
Stirling-cycle  compressors  introduce  undesirable  vibrations, 
and  thermoelectric  coolers  are  highly  inefficient  at  the  liquid- 
nitrogen  temperatures  of  interest)  and  to  cool  small,  special¬ 
ized  electronic  or  high-temperature-superconductor  devices 
(which  can  justify  the  expense  of  designing  an  optical 
cooler).  In  the  laboratory,  the  cooling  load  often  is  simply  the 
absorbing  sample  itself.  In  that  case,  the  cooling  energy  first 
reduces  the  temperature  of  the  material  (suspended  in  an 
optical-access  vacuum  chamber)  and  eventually  balances  the 
heat  load  on  the  sample  from  the  surroundings.  In  a  well- 
designed  system,  this  load  arises  primarily  from  ambient 
thermal  radiation,  implying  that  all  three  inputs  and  outputs 
in  Fig.  1  are  radiative.  Optical  heat  shields  (sometimes  called 
hot  mirrors)  can  be  used  to  further  reduce  this  heat  leak.  The 
sides  of  the  active  sample  are  coated  with  a  material  that 


thermal 
excitation 


7K 


3 

2 


laser 

pump 


hvf 

^  mean 
fluorescence 


transmits  the  fluorescence  but  reflects  the  ambient  blackbody 
light  (assuming  that  the  fluorescence  and  the  blackbody  ra¬ 
diation,  which  peaks  at  10  /rm  in  the  case  of  300  K  surround¬ 
ings,  are  spectrally  well  separated  from  one  another). 

The  current  record  temperature  drop  for  an  unshielded 
solid  sample  suspended  in  vacuum  starting  from  room 
temperature7  is  an  impressive  65  °C.  This  drop  was  achieved 
for  optical  pumping  of  the  rare-earth  ion  Yb3  +  doped  into  a 
heavy-metal  fluoride  glass  (ZBLAN).  The  advantages  of 
Yb3  +  over  other  species  are  threefold:  It  consists  of  only  two 
bands  of  energy  levels,  thus  avoiding  excited-state  absorp¬ 
tion  of  the  pump  or  fluorescence  radiation;  the  two  bands  are 
each  split  into  a  number  of  closely  spaced  levels,  thus  per¬ 
mitting  efficient  absorption  of  thermal  energy  from  the  host; 
and  the  energy  gap  between  the  two  bands  is  large  (corre¬ 
sponding  to  a  near-infrared  wavelength  of  1  /xm),  thereby 
minimizing  nonradiative  decay.8  Although  helpful,  these 
three  factors  are  not  required  for  successful  cooling  of  rare- 
earth  ions,  as  evidenced  by  recent  results9  for  Tm3  +  .  Laser 
cooling  of  gaseous  carbon  dioxide,10  organic  laser  dyes  in  an 
alcohol  solution,1112  and  gallium  arsenide  hetero¬ 
structures1314  also  have  been  reported. 

The  alert  reader  may  notice  that  these  systems  (rare-earth 
doped  solids,  C02,  solvated  organic  dyes,  and  semiconduc¬ 
tor  heterostructures)  also  make  efficient  lasers.  The  proper¬ 
ties  required  for  a  material  to  perform  well  as  an  optical 
cooler  (efficient  luminescence,  large  optical  cross  sections, 
minimal  competing  decay  channels)  also  are  properties  that 
promote  lasing.  In  fact,  one  need  simply  reverse  every  arrow 
in  Fig.  1  to  obtain  the  optically  pumped  laser  sketched  in  Fig. 
3.  The  optical  pump  might  now  be  a  flashlamp,  which,  like 
the  fluorescence  in  Fig.  1,  is  broadband,  incoherent,  undi¬ 
rected,  and  unpolarized.  Some  fraction  of  this  radiative  en¬ 
ergy  is  converted  into  laser  light  and  the  rest  appears  as 
waste  heat  to  be  removed  say  by  circulating  cooling  water. 

The  realization  that  good  lasers  can  also  make  good  opti¬ 
cal  coolers  suggests  other  possible  cooling  candidates.  Both 
Nd3  +  :YAG  and  ruby  have  been  explored  for  this  role,15’16 
although  for  technical  reasons  these  particular  materials  have 
not  lived  up  to  their  theoretical  cooling  potential.17 


III.  REVIEW  OF  RADIATION  THERMODYNAMICS 


Fig.  2.  A  possible  energy-level  scheme  of  a  cooling  material.  The  three 
states  could  represent  electronic  levels  of  an  atom,  vibrational  levels  of  a 
molecule,  or  levels  within  the  valence  and  conduction  bands  of  a  semicon¬ 
ductor. 


The  rate  at  which  entropy  is  carried  by  a  steady,  unpolar¬ 
ized  beam  of  light  (not  necessarily  collimated)  in  vacuum 
across  a  surface  A  is  given  by 
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S  =  2kc 


2 


[( 1  +«)ln(l  +n) 


(1) 


3 


A  J  Cl  J  Av 


—  n  In  n  ]  v2  d  v  cos  Odd  dA, 


as  discussed  in  the  Appendix,  where  c  is  the  speed  of  light, 
and  6  and  <f>  are  the  polar  and  azimuthal  angles,  respectively, 
at  which  a  photon  is  traveling  into  element  d£l 
=  sin  OdOdcf)  of  a  solid  angle  relative  to  the  normal  to  an 
element  dA  of  the  surface.  (Note  that  A  could  be  part  or  all 
of  an  arbitrary  cross  section  of  the  beam,  or  it  could  be  a 
portion  of  the  real  surface  of  a  material  emitting  or  absorbing 
the  light.)  The  radiation  is  distributed  over  a  set  of  optical 
modes  with  mean  occupation  number  n  that  depends  on  the 
photon  frequency  v,  the  direction  of  travel  6  and  <b  of  the 
photon  (within  a  range  of  solid  angle  O),  and  two  position 
coordinates  of  the  photon  on  the  surface  A  (affording  a  con¬ 
siderable  opportunity  for  confusion  if  spherical  coordinates 
are  used  to  describe  it).  The  frequency  integration  is  over  the 
relevant  spectral  bandwidth  An  of  the  light,  where  the  factor 
of  v 2  comes  from  the  density  of  states.  Finally  the  term  in 
square  brackets  in  Eq.  (1)  arises  from  Bose-Einstein  statis¬ 
tics.  It  is  crucial  to  realize  that  Eq.  (1)  is  valid  whether  or  not 
the  radiation  is  thermal,  that  is,  even  for  a  nonequilibrium 
photon  distribution  such  as  the  fluorescence  emission  of 
Fig.  2. 

Entropy  leads  to  the  concept  of  radiation  temperature.  Two 
definitions  of  this  quantity  have  been  proposed,18  called  the 
brightness  and  flux  temperatures.  Because  the  energy  of  a 
photon  is  hv  and  the  average  number  of  such  photons  in  a 
particular  beam  mode  is  n,  the  replacement  of  k  times  the 
quantity  in  square  brackets  in  Eq.  (1)  by  nhv  gives  the  rate  at 
which  energy  is  carried  by  the  beam  (that  is,  its  power), 

E  =  2hc~2  I  I  I  nv3dv cos  6 dH  dA.  (2) 

J  A  J  Cl  J  Av 


From  Eq.  (2)  the  spectral  radiance  (often  referred  to  as 
brightness,19  although  properly  speaking  that  is  a  photomet¬ 
ric  not  a  radiometric  term)  Lv=dE/(dA±dv  dEl),  where 
dA±  =  cos  OdA,  is  seen  to  be 

Lv=2nhv3/c2.  (3) 


X  (nm) 


Fig.  4.  The  spectral  radiance  of  a  5800  K  (solar  surface  temperature)  black- 
body  and  a  Gaussian  source  with  a  peak  at  \0=  1  /xm  representing  a  LED 
with  a  hemispherical  irradiance  of  1  W/mm2  and  a  25  nm  bandwidth  (w 
—  10  nm).  Within  this  bandwidth,  the  area  under  the  Planck  curve  is  equal 
to  the  area  under  the  Gaussian.  We  could  therefore  describe  this  LED  as 
being  as  bright  as  the  sun. 


Here  o-=2r75k4/(15c2/z3)  =  5.67X  1CT8  WnT2  K”4  is  the 
Stefan-Boltzmann  constant. 

Likewise  its  entropy  flux  density  Is=dS/dA  is18 


tBB 

Js 


2t rk  f°° 

c2  Jo 


h  v/kT B 

1  —  exp  (  —  hv/kTB) 


—  \n[exp(hv/  kTB)  —  1] 


X  v2dv=  —crTg. 


(6) 


The  thermodynamic  definition  of  temperature  is 

1  _  dS 
T=dE 


(7) 


If  we  apply  Eq.  (7)  to  the  cavity  radiation  escaping  from  the 
interior  of  a  furnace  of  volume  V,  we  find 


T= 


dr. 


dir 


AcrTBdTB 

AcrT2BdTB 


=  T, 


(8) 


Thus  an  experimental  measurement  of  the  spectral  radiance 
(using  an  angle-resolved  absolute-intensity  spectrometer)  di¬ 
rectly  gives  the  photon  occupation  numbers  n . 

The  brightness  temperature  TB  of  a  beam  of  radiation  is 
defined  as  the  temperature  of  a  blackbody  whose  spectral 
radiance  integrated  over  the  bandwidth  of  the  radiation  is 
equal  to  the  (integrated)  radiance  of  the  beam, 

nv3dv=\  -  , — ~dv,  (4) 

JAv  JA,exp(hv/kTB)-l 


In  agreement  with  Eq.  (4),  the  brightness  temperature  of 
blackbody  radiation  is  equal  to  the  absolute  temperature  of 
the  furnace  emitting  that  radiation.  Hence  the  irradiance  of 
blackbody  radiation  is  only  a  function  of  the  temperature  of 
its  source,  and  it  is  in  this  sense  that  one  describes  blackbody 
radiation  as  “pure”  heat  or  “low-grade”  energy. 

Consider  next  the  case  of  narrowband  radiation  centered  at 
peak  frequency  v0 ,  such  as  might  be  emitted  by  an  LED  or 
a  real  laser.  It  follows  from  Eq.  (4)  that  its  mean  brightness 
temperature  is 


according  to  Planck’s  law.  Further  insight  into  the  meaning 
of  the  brightness  temperature  can  be  obtained  by  examining 
two  special  cases  of  radiation.  First  consider  a  blackbody 
emitting  isotropically  into  a  hemisphere  (via  a  small  hole  in 
a  furnace  for  example).  Its  energy  flux  density  (or  irradiance) 

IE=dE/dA  is20 


tBB 

le 


2ttIi  f00  v3  4 

c2  Jo  exp(hv/kTB)  —  1  a^B' 


(5) 


f»  =  tln(i?i/,T)-  (9) 

where  n  =  c2L vl (2hv3Q)  according  to  Eq.  (3).  Here  Lv  is  the 
spectral  radiance  of  the  source  averaged  over  its  bandwidth 
Av. 

The  intuitive  meaning  of  Eq.  (9)  is  easier  to  grasp  from  a 
graph  such  as  Fig.  4.  Suppose  the  source  spectrum  has  a 
Gaussian  distribution 
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^\  =  ^peakexp 


(10) 


where  Lpeak  is  the  spectral  radiance  at  the  source’s  peak 
wavelength  k0=c/v0  and  w  is  the  spectral  width  corre¬ 
sponding  to  one  standard  deviation.  If  we  define  the  source’s 
bandwidth  as 


A\  = 


r. 


(ID 


where  g(A)^Lk/Lpeak  is  the  normalized  profile,  then  one 
can  show  that 

A\  =  w  \J2tt=2w ^2  In 2,  (12) 


which  is  the  full-width  at  half-maximum  of  the  source’s 
spectrum.  It  is  left  as  an  exercise  to  verify  that 


2 hc2\0 5 

Lpe^LBB(A0)  = - = - . 

exp  (hdkT  b\q)—  1 


(13) 


That  is,  to  an  excellent  approximation  one  can  determine  the 
brightness  temperature  of  a  narrowband  source  by  plotting 
its  spectral  radiance  and  then  finding  the  blackbody  curve 
which  passes  through  its  peak,  as  in  Fig.  4. 

For  a  very  bright  source,  Eq.  (9)  becomes  kT B=nhv0 
=  XqZ„/2.  For  example,  an  unpolarized  1  mW  red  helium- 
neon  laser  with  a  beam  area  of  1  mm2,  a  divergence  (half¬ 
angle  S)  of  0.5  mrad  corresponding  to  a  solid  angle  of  ap¬ 
proximately  ttS2  =  0.8  /jl sr,  and  a  bandwidth  of  1  GHz  has  a 
mean  brightness  temperature  on  the  order  of  2X1010  K. 19 
(The  detailed  variation  in  the  brightness  temperature  with 
frequency  and  angle  can  easily  be  deduced  for  simple  pro¬ 
files  such  as  Gaussians.21) 

If  the  narrowband  radiation  is  independent  of  the  angular 
directions  0  and  c[>  over  a  circular  cone  of  half-angle  S ,  then 
its  energy  flux  density  is 

I™=2Trhc~2nvoAv  sin2  <5,  (14) 


and  its  entropy  flux  density  is 

I^B=2TTkc~2[(  1  +n)ln(  1  +n)  —  n  InnJi'gAr'sin2  S. 

(15) 


The  analog  of  Eq.  (8)  becomes 

dI^B  hv0dn 
1  dI^B  k  ln(  1  +  1  ln)dn 


=  Tb 


(16) 


using  Eq.  (9)  in  the  last  step.  Therefore  the  brightness  tem¬ 
perature  remains  an  absolute  thermodynamic  temperature  for 
nonequilibrium  radiation. 

It  is  worth  noting  what  happens  if  either  the  bandwidth  Av 
or  divergence  S  of  the  source  collapses  to  zero.  If  its  energy 
flux  density  is  to  remain  finite,  then  Eq.  (14)  implies  that 
n— > oo.  In  this  limit,  Eq.  (15)  becomes 

I^B=2Trkc~2\(  1  +n)lnn  —  n  In if]  A  v  sin2  S 


_L,nb!^ 

I  *  F  — 

n  J'o  w 


0. 


(17) 


That  is,  the  entropy  carried  by  monochromatic  or  unidirec¬ 
tional  radiation  is  zero,  so  that  one  can  characterize  an  ideal 
laser  beam  as  “pure”  work  or  “high-grade”  energy.  Note 
from  Eq.  (16),  however,  that  JT  dls  remains  equal  to  the 


irradiance  (just  as  it  did  for  blackbody  radiation)  and  is  not 
zero,  as  one  might  have  expected  if  this  integral  defined  a 
heat  flux  density.  Entropy  and  heat  are  not  directly  related  for 
nonequilibrium  distributions. 

To  analyze  the  efficiency  of  optical  devices,  it  is  useful  to 
introduce  the  flux  temperature, 

TF=y.  (18) 

lS 

Equations  (14)  and  (15)  imply  for  a  narrowband  source  that 
h  vn  n 

T  = _  1191 

F  k  ( 1  +n)ln(  1  +n)  —  n  In  u  ' 


If  the  source  is  no  more  than  moderately  bright,  then  n<  1, 
and  Eq.  (19)  becomes 

h  Vn  n  li  Vn  /  k 

j  3 _ _ _ = _ _ _  /'orq 

F  k  (\){n)  —  n\nn  l+ln(l/n)' 

Equations  (9)  and  (20)  imply  that  TbITf—>  1  as  n— >  0,  so 
that  the  flux  and  brightness  temperatures  are  equal  for  dim, 
narrowband  radiation. 

In  contrast,  TF=0.15TB  for  blackbody  radiation  according 
to  Eqs.  (5)  and  (6).  At  first  glance  this  result  is  surprising, 
because  it  implies  that  if  a  large  body  (of  surface  area  A)  at 
temperature  T  radiates  away  a  small  amount  of  heat  Q  into 
free  space,  then  the  light  carries  away  entropy  .S' =  4  (1/3  7’ 
according  to  Eq.  (18),  rather  than  QIT.  This  entropy  trans¬ 
port  appears  to  disagree  with  the  first  equality  in  Eq.  (8), 
which  can  be  rewritten  in  this  context  as  S  =  JdQIT.  But  in 
fact  Eq.  (8)  does  not  imply  that  S  =  fdQ/T=Q/T,  because 
Q  is  only  a  function  of  temperature  (specifically  dQ 
=  4crAT2dT),  whereas  Q  is  implicitly  being  interpreted  as  a 
function  of  time  at  the  (nearly)  constant  temperature  of  the 
emitter.  The  lesson  is  to  exercise  caution  when  relating  the 
fluxes  to  the  total  energy  and  entropy.  The  light  carries  away 
entropy  4(1/373  while  the  entropy  of  the  body  decreases  by 
QIT ,  and  thus  the  total  entropy  change  of  the  universe  is 
positive  in  this  irreversible  emission  process. 

On  the  other  hand,  we  expect  entropy  to  balance  if  a  hot 
body  at  temperature  7’+A7’  is  radiatively  coupled  to  sur¬ 
roundings  at  temperature  T  that  are  only  infinitesmally 
smaller  in  temperature  by  A  73  In  this  case,  the  net  rate  of 
radiative  energy  transport  from  the  hot  body  to  the  surround¬ 
ings  is 


Q  =  o-A\(T+AT)4-T4]  =  4o-AT3  AT,  (21) 

while  the  net  rate  of  entropy  transport  is 
4 

S=  —crA\(T  +  AT)3  —  T3]=4 crAT2  AT.  (22) 


This  time  S  =  Q/T,  as  one  would  predict  from  the  quasi¬ 
equilibrium  of  the  bodies  with  the  radiation. 

Returning  to  our  discussion  of  Eq.  (19),  for  a  very  bright 
source  («8>1)  we  obtain 

T  J1V o  n  _,lvo  n 

F  k  (1  +n)lnn  —  n  Inn  k  In n' 

For  ideal  laser  radiation,  n— >° °  and  thus  T F  becomes  infinite, 
which  is  consistent  with  the  fact  that  it  carries  zero  entropy  at 
a  finite  irradiance.  One  (but  not  the  only)  way  to  get  zero 
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entropy  transport  is  if  the  radiation  occupies  a  single  optical 
mode,  corresponding  to  a  monochromatic  plane  wave,  so 
that  the  multiplicity  of  its  macrostate  is  unity.  Also  note  that 

Eqs.  (9)  and  (23)  imply  that  TB/TF=\nn  for  a  bright  nar¬ 
rowband  source.  That  is,  the  brightness  temperature  diverges 
even  more  rapidly  than  the  flux  temperature  as  n— 

So  far  it  has  been  implicitly  assumed  that  the  fluorescent 
emission  in  Fig.  1  is  into  free  space.  However  if  the  cooling 
sample  is  bathed  in  thermal  radiation  from  the  surroundings 
at  ambient  temperature  T A ,  then  one  needs  to  account  both 
for  the  occupation  number  n  of  fluorescence  photons  and  the 
distribution  nA  of  ambient  photons.22  The  net  energy  and 
entropy  fluxes  stem  from  the  difference  between  the  radia¬ 
tion  leaving  the  body  (fluorescence  plus  thermal  radiation 
which  is  reflected  or  transmitted,  assuming  for  simplicity  that 
the  sample  has  negligible  thermal  emissivity)  and  that  inci¬ 
dent  on  the  body  (thermal  radiation).2’  Hence 


E„et=2hc~ 


!J  (n  +  n 


A)v3dv  cos  0dfl  dA 


—  2hc  2J  nAv3dv  cos  6  dfl  dA 
=  2 hc~2  J  nv3dv cos  0 dTl  dA, 


(24) 


so  that  one  recovers  Eq.  (2).  In  contrast 


Snet=2kc  2J  [(1  +  n  +«A)ln(l  +n  +  nA) 

—  (n  +  nA)\n(n  +  nA)  —  (l  +nA)ln(l  +nA) 

+  nA\nnA]v2dv cos  0 d£l  dA,  (25) 


which  does  not  simplify  to  Eq.  (1).  That  is,  the  entropy  car¬ 
ried  by  the  net  radiative  emission  from  the  refrigerator  is  not 
independent  of  the  ambient  temperature  in  general.  However 
for  practical  optical  refrigerators,  n>nA,  and  Eq.  (1)  is  an 
excellent  approximation  to  the  net  entropy  flux.  It  is  only  in 
the  limit  of  vanishing  fluorescence  radiance  that  the  ambient 
radiation  matters.  Specifically  Eq.  (25)  can  then  be  approxi¬ 
mated  as 


S nei=  2kc  2J  [(1  +/r  +  nA)ln(l  +nA)~  (n  +  nA)lnnA 
—  ( 1  +nA) ln(  1  +nA)  +  nA  In  nA]  v2dv  cos  0  dTl  dA 
=  2 kc  2  J  n  ln(  1  +  \lnA)v2dv  cos  0dfl  dA 
=  EnJTA  (26) 


after  substituting  the  Planck  expression  for  nA  in  the  last 
step.  Consequently  the  lower  limit  of  the  fluorescence  flux 
temperature  Enet/Snet  is  the  ambient  temperature  T A  as  the 
pump  power  in  Fig.  1  is  reduced  toward  zero.  Even  if  the 
laser  pump  is  shut  off,  there  remains  some  weak  fluorescence 
stimulated  by  the  absorption  of  ambient  thermal  radiation. 


IV.  LIMITING  EFFICIENCIES  OF  OPTICAL 
CONVERTERS 


The  stage  is  set  for  a  thermodynamic  analysis  of  the  effi¬ 
ciency  with  which  the  optical  refrigerator  in  Fig.  1  converts 
heat  into  light.  According  to  the  first  law  (energy  conserva¬ 
tion), 

Eovx=EmAQ.  (27) 

The  cooling  coefficient  of  performance  is  defined  in  the 
usual  way  for  a  refrigerator  as 

(28) 

£in 


The  maximum  value  of  77  is  the  Carnot  limit,  rjc,  and  is 
determined  by  the  second  law  of  thermodynamics.  The  en¬ 
tropy  carried  away  by  the  fluorescence  cannot  be  less  than 
the  sum  of  the  entropy  withdrawn  from  the  cooling  sample 
and  the  entropy  transported  in  by  the  pump  laser. 


^out  ^  Q  Em 

~y;^t  v 


(29) 


where  T  is  the  steady-state  operating  temperature  of  the  re¬ 
frigerator,  and  T f  and  Tp  are  the  flux  temperatures  of  the 
fluorescence  and  pump  radiation,  respectively.  The  reversible 
Carnot  limit  is  obtained  by  choosing  the  equality  sign  in  Eq. 
(29).  If  one  substitutes  Eqs.  (27)  and  (28)  into  Eq.  (29),  we 
obtain 


T- AT 

Vc=^r^F,  (30) 

Lf  1 

where  AT=TT f/T p  .  Equation  (30)  can  be  interpreted  as  the 
Carnot  coefficient  of  performance  of  the  engine-driven  re¬ 
frigerator  depicted  in  Fig.  5.  The  engine  extracts  heat  Em 
from  a  hot  reservoir  at  temperature  Tp  and  rejects  heat  fne 
to  a  cold  reservoir  at  temperature  T f.  It  thus  produces  work 
W  with  a  Carnot  efficiency  of 


^7engine 


(31) 


This  work  then  drives  a  refrigerator  to  draw  heat  Q  out  of  a 
cold  reservoir  at  temperature  T  and  to  dump  waste  heat 
£^”jge  into  a  hot  reservoir  at  temperature  Tf,  with  a  Carnot 
coefficient  of  performance  given  by 


^7  fridge 


(32) 


Note  that  as  the  pump  power  Em  is  reduced,  T ,  decreases — 
with  a  lower  limit  of  TA^T  according  to  Eq.  (26) — and 
hence  r]fAAge  increases,  but  at  the  expense  of  an  overall  de¬ 
creased  cooling  power  rjcEin.  The  product  of  Eqs.  (31)  and 
(32)  equals  Eq.  (30),  and  the  net  input  and  output  energy 
fluxes  of  Fig.  5  reproduce  those  of  Fig.  1  if  one  identifies 
Eout=ElTe+Efnage_  scheme  prompts  the  realization 
that  the  engine  is  required  only  because  the  pump  radiation  is 
not  “pure  work”  for  a  real  laser.  If  the  pump  laser  were  ideal, 
then  Tp— > 00  so  that  A  T=  0  and  Eq.  (30)  would  directly  re¬ 
duce  to  Eq.  (32).  In  contrast,  if  one  were  to  try  to  pump  an 
optical  cooler  with  fluorescence  radiation  (say  by  recycling 
the  light),  then  Tp—>Tj  and  Eq.  (30)  indicates  that  its  coef- 
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he 


1 


(35) 


Fig.  5.  A  heat  engine  operating  between  thermal  reservoirs  with  a  high 
temperature  of  Tp  and  a  cold  temperature  of  Tf  which  is  moderately  larger 
than  ambient,  whose  output  work  W  is  used  to  drive  a  refrigerator,  which 
cools  a  load  at  temperature  T  somewhat  less  than  ambient  and  rejects  heat 
into  the  reservoir  at  intermediate  temperature  T * . 


ficient  of  performance  would  fall  to  zero.  Thermodynamics 
thus  explains  why  optical  cooling  requires  a  laser  pump  for 
efficient  operation:  The  loss  in  entropy  of  the  cooling  load 
must,  at  minimum,  be  compensated  by  the  gain  in  entropy  of 
the  radiation  as  it  is  converted  from  the  input  pump  to  the 
output  fluorescence.  Therefore,  it  is  not  surprising  that 
Landau24  dismissed  the  practicality  of  fluorescent  cooling  in 
1946. 

Consider  an  example  using  actual  values22  relevant  to  la¬ 
ser  cooling  of  Yb3  +  :ZBLAN.  The  laser  pump  is  narrowband 
and  very  bright,  so  that  Eq.  (23)  is  applicable. 


lie  np 
k\p  In  np  ’ 


where 


(33) 


p  2 he  TTR~AvpTrSp 


(34) 


from  Eq.  (3).  Suppose  that  the  pump  laser  has  a  power  of 
Ep^40  W,  a  beam  radius  of  Rp  =  0.5  mm,  a  bandwidth  of 
A  1^  =  40  GHz,  a  divergence  of  Sp=  1  mrad,  and  a  wave¬ 
length  of  \p  =  1 030  nm.  Then  1 0<J,  which  justifies  the 
use  of  Eq.  (23),  and  Tp^7X  1011  K.  In  contrast,  the  fluores¬ 
cence  is  only  moderately  bright,  so  that  Eq.  (20)  is  called  for. 


Tf  k\0fl+]n(l/nf)’ 

with 

_  =  XQ/ 

2  h  c  A  sampie A  Vy 7T 

The  fluorescence  spectrum  has  a  peak  wavelength  of  \0y 
=  975  nm  and  a  full-width  at  half-maximum  of  about  A\f 
=  35  nmsc  Xq/Avj  at  room  temperature.  The  fluores¬ 
cence  is  assumed  to  be  emitted  homogeneously  and  hemi- 
spherically  from  the  surface  of  the  cooling  sample  (of  area 
A  sample)  >  which  is  taken  to  be  a  cylinder  whose  height  and 
diameter  are  each  3.0  cm,  with  a  total  power  of  Ef=  40  W. 
[Note  that  Ep  =  Eft(\  +  rj)  according  to  Eqs.  (27)  and  (28), 
thus  explaining  the  approximation  in  the  specification  of  the 
pump  power.]  Hence  Hy=6.4XlO~4  and  Zy=1760  K.  (An 
exact  calculation  using  a  measured  fluorescence  spectrum17 
gives  a  similar  value  of  1530  K,  thereby  justifying  the  ap¬ 
proximations  used  here.)  Clearly  Tp  is  so  much  larger  than 
T f  that  Eq.  (32)  can  be  taken  to  be  the  Carnot  coefficient  of 
performance  of  this  optical  refrigerator.  Consequently  r]c 
=  20%  at  room  temperature,  and  it  diminishes  approximately 
linearly  to  zero  as  71— >0. 

In  contrast,  the  actual  cooling  coefficient  of  performance 
is 


factual  =(vf-Vp)/vp=(\p-Af)/Af,  (37) 

according  to  the  discussion  in  Sec.  II,  assuming  that  the 
pump  beam  is  entirely  absorbed  by  the  sample.  In  the  case8 

of  Yb3  +  :  ZBLAN,  \f=  995  nm  so  that  factual =  3 -5%  at  the 
pump  wavelength  of  1030  nm  used  in  Eqs.  (33)  and  (34). 
Although  it  might  appear  that  one  could  increase  the  cooling 
performance  by  tuning  the  pump  laser  to  longer  wavelengths, 
in  practice  the  measured  value  of  rj  is  found  to  roll  off  be¬ 
cause  (heat  producing)  trace  impurity  absorption  begins  to 
dominate  the  rapidly  decreasing  ytterbium  absorption.  Fur¬ 
thermore,  the  cooling  coefficient  decreases  faster  than  lin¬ 
early  as  the  operating  temperature  is  decreased,22  because  the 
long-wavelength  absorption  coefficient  decreases  with  de¬ 
creasing  temperature  (so  that  \p  has  to  be  reduced),  while 
simultaneously  the  average  fluorescence  wavelength  in¬ 
creases  (due  to  diminishing  phonon  absorption  within  the 
ground  and  excited  bands  of  energy  levels). 

An  alternative  to  pumping  the  refrigerator  in  Fig.  1  opti¬ 
cally  is  to  pump  a  semiconductor  diode  electrically.25  The 
key  to  this  electroluminescent  cooling  is  to  choose  an  applied 
forward  bias  V  which  is  sufficiently  smaller  than  the  average 
bandgap  emission  energy  h  V  divided  by  the  electron  charge 
e .  As  in  the  case  of  laser  cooling,  one  factor  that  determines 
how  much  smaller  the  pump  energy  e  V  must  be  than  h  V  to 
achieve  net  cooling  is  the  fluorescence  quantum  efficiency. 
In  particular.  Auger  processes  and  surface  recombination 
limit  this  efficiency.  If  we  ignore  such  limitations,  the  Carnot 
coefficient  of  performance  is  again  given  by  Eq.  (32),  be¬ 
cause  the  electrical  pump  in  this  case  is  a  source  of  “pure” 
work.  In  practice,  however,  the  required  electrical  current 
creates  Joule  heating  of  the  diode,  which  is  the  primary  rea¬ 
son  that  optical  pumping  is  preferred.13  Unfortunately,  the 
large  refractive  index  of  semiconductors  (e.g.,  3.6  for  GaAs) 
leads  to  trapping  of  the  fluorescence  photons,  which  de- 
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creases  the  external  fluorescence  quantum  efficiciency,  as  the 
photons  rattle  around  inside  the  sample  and  are  thus  more 
likely  to  encounter  a  nonradiative  decay  channel.  Berdahl26 
has  suggested  that  this  total  internal  reflection  could  be  frus¬ 
trated  by  bringing  an  external  absorber  to  within  a  few  mi¬ 
crons  of  the  semiconductor  surface.  Another  suggestion  by 
Berdahl  is  to  reverse  bias  the  diode  so  that  it  cools  a  neigh¬ 
boring  sample  by  absorbing  its  thermal  emission,  a  phenom¬ 
enon  known  as  negative  luminescence. 

The  preceding  thermodynamic  analysis  of  photolumines- 
cent  cooling  can  be  similarly  applied  to  the  optically  pumped 
laser  in  Fig.  3.  We  define  its  efficiency  as 


(38) 


An  analysis  similar  to  the  derivation  of  Eq.  (30)  can  be  used 
to  deduce  a  Carnot  limit  of 


Vc= 


T  —  T 

1  p  1 

Tp—  AT’ 


(39) 


where  A  T=  TTp  /  Tt .  In  this  case,  T  is  the  steady-state  oper¬ 
ating  temperature  of  the  laser  medium,  and  T /  and  Tp  are  the 
flux  temperatures  of  the  output  laser  and  input  pump  radia¬ 
tion,  respectively.  As  one  might  have  guessed,  rjc  increases 
as  the  ratio  TpITi  of  these  two  flux  temperatures  increases.27 
In  fact,  the  thermodynamic  limit  on  the  efficiency  of  the  laser 
becomes  100%  as  this  ratio  approaches  unity,  which  is  a  big 
advantage  of  diode-pumped  laser  systems.  (In  practice,  how¬ 
ever,  Tp  will  always  be  less  than  T ,  because  a  set  of  mutually 
incoherent  diode  laser  beams  is  essentially  being  combined 
into  a  single,  high  quality  output  beam.) 

The  idea  of  pumping  one  laser  with  another  suggests  an 
intriguing  concept.  One  can  connect  the  refrigerator  of  Fig.  1 
to  the  laser  of  Fig.  3,  and  simultaneously  pump  both  in  such 
a  manner  that  all  of  the  waste  heat  generated  by  the  latter  is 
removed  by  the  former.  Remarkably,  this  scheme  can  be 
implemented  using  the  same  active  material  for  both 
devices,28  by  choosing  the  pump  frequency  vp  to  be  interme¬ 
diate  between  the  mean  spontaneous  emission  frequency  Tj 
and  the  stimulated  emission  frequency  .  In  effect,  one  is 
balancing  the  anti-Stokes  cooling  from  the  refrigerator 
against  the  Stokes  heating  from  the  laser,  resulting  in  what  is 
known  as  athermal  or  radiation-balanced  lasing.  A  simulta¬ 
neous  optical  and  thermodynamic  analysis  of  a  model  ather¬ 
mal  laser  based  on  an  ytterbium-doped  crystalline  system 
can  be  found  in  Ref.  29. 


V.  CONCLUDING  REMARKS 

Laser  cooling  is  now  a  well  established  laboratory  tech¬ 
nique  with  myriad  practical  applications.  For  atom  cooling, 
the  concept  frequently  is  explained  from  the  point  of  view  of 
momentum,  in  which  large  numbers  of  photons  are  fired  like 
high-speed  ping-pong  balls  against  an  oncoming  bowling- 
ball-like  atom  to  translationally  slow  it  down.30  Another  ex¬ 
planation,  appropriate  both  to  the  idea  of  the  anti-Stokes 
cooling  of  solids  and  Doppler  cooling  of  atoms,  is  to  use 
energy  conservation  to  balance  the  heat  loss  from  the  target 
against  the  energy  gain  of  the  radiation.  As  good  as  these 
explanations  are,  neither  is  able  to  explain  why  a  laser  pump 
is  required  rather  than  say  a  collimated  beam  of  narrowband 
light  from  a  high-intensity  lamp.  The  simplest  explanation  of 
this  last  issue  invokes  entropy  as  discussed  in  this  paper. 


In  particular,  a  statistical  equation  for  radiation  entropy 
can  be  deduced  from  the  Shannon  information- theoretical 
definition.  The  advantage  of  this  approach  over  the  tradi¬ 
tional  radiative  heat  transfer  viewpoint  is  that  it  encompasses 
nonequilibrium  photon  distributions.  Once  the  entropy  is  de¬ 
fined,  two  effective  temperatures  of  radiation  can  be  intro¬ 
duced  to  parametrize  it.  A  Camot  limit  on  the  performance  of 
an  optical  engine  or  refrigerator  is  then  easily  derived  from 
the  first  and  second  laws  of  thermodynamics.  As  usual,  this 
limit  would  only  be  attainable  for  reversible  operation  of  the 
device.  The  calculations  in  this  article  indicate,  for  example, 
that  an  ytterbium  based  solid-state  laser  cooler  has  a  thermo¬ 
dynamic  limit  on  its  coefficient  of  performance  of  20%  for  a 
typical  set  of  operating  conditions.  However,  experimental 
measurements  to  date  have  indicated  that  at  best  3.5%  of  the 
pump  power  has  actually  been  converted  into  cooling  using 
this  material.  This  discrepancy  is  indicative  of  the  inefficien¬ 
cies  resulting  from  impurities  in  the  samples  and  irrevers¬ 
ibilities  in  the  fluorescent  emission.  On  the  positive  side,  this 
shortfall  implies  that  the  ultimate  limits  on  optical  cooling  of 
bulk  matter  have  not  been  approached  and  that  plenty  of 
room  remains  for  new  scientific  and  engineering  develop¬ 
ments  in  this  emerging  field  of  research. 1 
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APPENDIX:  THE  ENTROPY  OF  RADIATION 

A  recent  derivation  of  Eq.  (1)  can  be  found  in  Ref.  21,  but 
it  is  beyond  the  understanding  of  undergraduate  students.  A 
more  accessible  treatment  is  in  Ref.  18  and  is  outlined  here. 
Filling  in  the  details  in  the  derivation  makes  for  a  reasonable 
student  exercise.  We  begin  from  the  Shannon  definition  of 
the  entropy, 

S=~k  X  P state Itl P state  ’  (Al) 

states 

where  the  summation  is  over  all  possible  states  of  the  sys¬ 
tem,  each  of  which  has  probability  P  state .  In  the  present  con¬ 
text,  a  particular  state  is  defined  by  the  occupation  of 
photons  in  optical  mode  1,  N2  photons  in  optical  mode  2, 
and  so  on.  Now  assume  that  the  probability  p,(A,)  of  finding 
N i  photons  in  mode  i  is  independent  of  the  probability 
Pj(Nj)  of  finding  Nj  photons  in  mode  j  whenever  i¥=j. 
Then 

00 

^ state ~  II  pm,  (A2) 

/=1 

where 

00 

2  Pi(Nj)=l  for  any  i.  (A3) 

Nj  =  0 

Straightforward  algebra  leads  to 

oo  oo 

s=-k 2  2  pm^pm-  (A4> 

i=  1  Nt  =  0 

Next  assume  that  the  probability  of  finding  one  additional 
photon  in  any  mode  is  independent  of  the  number  already 
occupying  that  mode.  This  assumption  implies  that 
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Pi(Ni )««?',  (A5) 

where  0^qt<  1.  If  we  normalize  Pi(Nj)  in  Eq.  (A5)  accord¬ 
ing  to  Eq.  (A3),  we  obtain 

Pi{Ni)  =  {\-qi)qNii .  (A6) 

But  the  definition  of  the  mean  occupation  number  is 

00 

nt=  2  -A iPi(N ,).  (A7) 

Nj  =  0 

We  substitute  Eq.  (A6)  into  Eq.  (A7)  and  evaluate  the  sum  to 
obtain 


<37 


1  +11;' 


(A8) 


Equation  (A8)  is  consistent  with  the  fact  that  q ,  can  be  in¬ 
terpreted  from  Eq.  (A6)  as  either  the  relative  probability  of 
there  being  one  more  photon  in  mode  i,  or  the  probability  of 
finding  a  nonzero  number  of  photons  in  that  mode.  (Consider 
for  example  the  cases  where  «,  =  0,  1,  or  °°.)  We  next  substi¬ 
tute  Eq.  (A8)  into  Eq.  (A6)  to  obtain 


p!{Ni)  = 


(l+«i) 


(A9) 


and  then  substitute  Eq.  (A9)  into  (A4)  to  finally  find 


00 


S  =  k^  [(1  +n,-) ln(l  -\-ni)  —  ni  In  nz-]. 

i=  l 


(A10) 


The  density  of  photon  states20  per  polarization  mode  in  a 
frequency  interval  dv  and  element  of  solid  angle  dfl  is 
c~3 4 5 6 7 8v2dvdil.  We  multiply  this  result  by  two  because  of  the 
two  independent  transverse  directions  of  the  polarization  of 
light,  and  by  the  speed  c  to  obtain  the  entropy  per  unit  time 
rather  than  per  unit  distance.  We  also  multiply  by  the  element 
of  surface  area  dA,  as  well  as  by  cos  0  to  project  this  area 
element  perpendicularly  to  the  direction  of  photon  propaga¬ 
tion.  We  thereby  convert  the  sum  in  Eq.  (A  10)  into  the  inte¬ 
gral  given  by  Eq.  (1). 
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